Epigenetic disruptions in the offspring hypothalamus in response to maternal infection.

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

Histone deficiency and hypoacetylation in the aging retinal pigment epithelium.

Histones serve as a major carrier of epigenetic information in the form of post-translational modifications which are vital for controlling gene expression, maintaining cell identity, and ensuring proper cellular function. Loss of histones in the aging genome can drastically impact the epigenetic landscape of the cell leading to altered chromatin structure and changes in gene expression profiles. In this study, we investigated the impact of age-related changes on histone levels and histone acetylation in the retinal pigment epithelium (RPE) and retina of mice. We observed a global reduction of histones H1, H2A, H2B, H3, and H4 in aged RPE/choroid but not in the neural retina. Transcriptomic analyses revealed significant downregulation of histones in aged RPE/choroid including crucial elements of the histone locus body (HLB) complex involved in histone pre-mRNA processing. Knockdown of HINFP, a key HLB component, in human RPE cells induced histone loss, senescence, and the upregulation of senescence-associated secretory phenotype (SASP) markers. Replicative senescence and chronological aging in human RPE cells similarly resulted in progressive histone loss and acquisition of the SASP. Immunostaining of human retina sections revealed histone loss in RPE with age. Acetyl-histone profiling in aged mouse RPE/choroid revealed a specific molecular signature with loss of global acetyl-histone levels, including H3K14ac, H3K56ac, and H4K16ac marks. These findings strongly demonstrate histone loss as a unique feature of RPE aging and provide critical insights into the potential mechanisms linking histone dynamics, cellular senescence, and aging.

Use of a Candida albicans SC5314 PacBio HiFi reads dataset to close gaps in the reference genome assembly, reveal a subtelomeric gene family, and produce accurate phased allelic sequences.

Candida albicans SC5314 is the most-often used strain for molecular manipulation of the species. The SC5314 reference genome sequence is the result of considerable effort from many scientists and has advanced research into fungal biology and pathogenesis. Although the resource is highly developed and presented in a phased diploid format, the sequence includes gaps and does not extend to the telomeres on its eight chromosome pairs. Accurate SC5314 genome assembly is complicated by the presence of extensive repeated sequences and considerable allelic length variation at some loci. Advances in genome sequencing technology provide the tools to obtain highly accurate long-read data that span even the most-difficult-to-assemble genome regions. Here, we describe derivation of a PacBio HiFi data set and creation of a collapsed haploid telomere-to-telomere assembly of the SC5314 genome (ASM3268872v1) that revealed previously unknown features of the strain. ASM3268872v1 subtelomeric distances were up to 19 kb larger than in the reference genome and revealed a family of highly conserved DNA helicase-encoding genes at 10 of the 16 chromosome ends. We also describe alignments of individual HiFi reads to deduce accurate diploid sequences for the most notoriously difficult-to-assemble C. albicans genes: the agglutinin-like sequence (ALS) gene family. We provide a tutorial that demonstrates how the HiFi reads can be visualized to explore any region of interest. Availability of the HiFi reads data set and the ASM3268872v1 comparative guide assembly will streamline research efforts because accurate diploid sequences can be derived using simple in silico methods rather than time-consuming laboratory-bench approaches.

An improved Lodderomyces elongisporus NRRL YB-4239 genome assembly substantiated by its electrophoretic karyotype.

Pacific Biosciences long-read sequencing was used to improve the genome assembly for Lodderomyces elongisporus strain NRRL YB-4239 (ATCC 11503). The new assembly included eight chromosomes that were substantiated by the electrophoretic karyotype. The nuclear genome was 16.1 Mb (37.2% GC) with 5,740 genes predicted.

Genome sequence of Emydomyces testavorans type strain 16-2883 (ATCC TSD-145), an onygenalean fungus associated with freshwater aquatic turtle shell lesions.

Emydomyces testavorans is an onygenalean keratinophilic fungus associated with shell and skin lesions in freshwater aquatic turtles. The genome sequence presented here includes five contigs (ranging in size from 2.8 to 9.8 Mb; 31.8 Mb total; 40% GC) and a 92.2-kb mitochondrial genome. The nuclear genome predicted 7,550 genes.

Complete genome sequence of the archetype bile acid 7α-dehydroxylating bacterium, Clostridium scindens VPI12708, isolated from human feces, circa 1980.

Clostridium scindens strain VPI12708 serves as model organism to study bile acid 7α-dehydroxylating pathways. The closed circular genome of C. scindens VPI12708 was obtained by PacBio sequencing. The genome is composed of 3,983,052 bp, with 47.59% G + C, and 3,707 coding DNA sequences are predicted.

A Chromosome-Level Genome Assembly for Yamadazyma tenuis ATCC 10573.

Pacific Biosciences (PacBio) long-read sequencing was used to generate a chromosome-level genome assembly for Yamadazyma tenuis strain ATCC 10573. The assembly featured 7 chromosomes that matched the electrophoretic karyotype and a 26.5-kb circular mitochondrial genome. The nuclear genome was 10.8 Mb, with a GC content of 43%, and 5,340 predicted genes.

First Draft Genome Assembly of Root-Lesion Nematode Pratylenchus scribneri Generated Using Long-Read Sequencing.

Root-lesion nematodes (genus Pratylenchus) belong to a diverse group of plant-parasitic nematodes (PPN) with a worldwide distribution. Despite being an economically important PPN group of more than 100 species, genome information related to Pratylenchus genus is scarcely available. Here, we report the draft genome assembly of Pratylenchus scribneri generated on the PacBio Sequel IIe System using the ultra-low DNA input HiFi sequencing workflow. The final assembly created using 500 nematodes consisted of 276 decontaminated contigs, with an average contig N50 of 1.72 Mb and an assembled draft genome size of 227.24 Mb consisting of 51,146 predicted protein sequences. The benchmarking universal single-copy ortholog (BUSCO) analysis with 3131 nematode BUSCO groups indicated that 65.4% of the BUSCOs were complete, whereas 24.0%, 41.4%, and 1.8% were single-copy, duplicated, and fragmented, respectively, and 32.8% were missing. The outputs from GenomeScope2 and Smudgeplots converged towards a diploid genome for P. scribneri. The data provided here will facilitate future studies on host plant-nematode interactions and crop protection at the molecular level.

Draft genome assemblies of the avian louse Brueelia nebulosa and its associates using long-read sequencing from an individual specimen.

Sequencing high molecular weight (HMW) DNA with long-read and linked-read technologies has promoted a major increase in more complete genome sequences for nonmodel organisms. Sequencing approaches that rely on HMW DNA have been limited to larger organisms or pools of multiple individuals, but recent advances have allowed for sequencing from individuals of small-bodied organisms. Here, we use HMW DNA sequencing with PacBio long reads and TELL-Seq linked reads to assemble and annotate the genome from a single individual feather louse (Brueelia nebulosa) from a European Starling (Sturnus vulgaris). We assembled a genome with a relatively high scaffold N50 (637 kb) and with BUSCO scores (96.1%) comparable to louse genomes assembled from pooled individuals. We annotated a number of genes (10,938) similar to the human louse (Pediculus humanus) genome. Additionally, calling phased variants revealed that the Brueelia genome is more heterozygous (∼1%) then expected for a highly obligate and dispersal-limited parasite. We also assembled and annotated the mitochondrial genome and primary endosymbiont (Sodalis) genome from the individual louse, which showed evidence for heteroplasmy in the mitogenome and a reduced genome size in the endosymbiont compared to its free-living relative. Our study is a valuable demonstration of the capability to obtain high-quality genomes from individual small, nonmodel organisms. Applying this approach to other organisms could greatly increase our understanding of the diversity and evolution of individual genomes.

RNA editing of microtubule-associated protein tau circular RNAs promotes their translation and tau tangle formation.

Aggregation of the microtubule-associated protein tau characterizes tauopathies, including Alzheimer's disease and frontotemporal lobar degeneration (FTLD-Tau). Gene expression regulation of tau is complex and incompletely understood. Here we report that the human tau gene (MAPT) generates two circular RNAs (circRNAs) through backsplicing of exon 12 to either exon 7 (12→7 circRNA) or exon 10 (12→10 circRNA). Both circRNAs lack stop codons. The 12→7 circRNA contains one start codon and is translated in a rolling circle, generating a protein consisting of multimers of the microtubule-binding repeats R1-R4. For the 12→10 circRNA, a start codon can be introduced by two FTLD-Tau mutations, generating a protein consisting of multimers of the microtubule-binding repeats R2-R4, suggesting that mutations causing FTLD may act in part through tau circRNAs. Adenosine to inosine RNA editing dramatically increases translation of circRNAs and, in the 12→10 circRNA, RNA editing generates a translational start codon by changing AUA to AUI. Circular tau proteins self-aggregate and promote aggregation of linear tau proteins. Our data indicate that adenosine to inosine RNA editing initiates translation of human circular tau RNAs, which may contribute to tauopathies.

Rewritable two-dimensional DNA-based data storage with machine learning reconstruction.

DNA-based data storage platforms traditionally encode information only in the nucleotide sequence of the molecule. Here we report on a two-dimensional molecular data storage system that records information in both the sequence and the backbone structure of DNA and performs nontrivial joint data encoding, decoding and processing. Our 2DDNA method efficiently stores images in synthetic DNA and embeds pertinent metadata as nicks in the DNA backbone. To avoid costly worst-case redundancy for correcting sequencing/rewriting errors and to mitigate issues associated with mismatched decoding parameters, we develop machine learning techniques for automatic discoloration detection and image inpainting. The 2DDNA platform is experimentally tested by reconstructing a library of images with undetectable or small visual degradation after readout processing, and by erasing and rewriting copyright metadata encoded in nicks. Our results demonstrate that DNA can serve both as a write-once and rewritable memory for heterogenous data and that data can be erased in a permanent, privacy-preserving manner. Moreover, the storage system can be made robust to degrading channel qualities while avoiding global error-correction redundancy.

Rat Spinal Cord Injury Associated with Spasticity Leads to Widespread Changes in the Regulation of Retained Introns.

To determine molecular changes that correlate with long-term physiological changes after spinal cord injury associated with spasticity, we used a complete transection model with an injury at sacral spinal level S2, wherein tail spasms develop in rats weeks to months post-injury. Using Illumina and nanopore sequencing, we found that from 12,266 expressed genes roughly 11% (1,342) change expression levels in the rats with spasticity. The transcription factor PU.1 (Spi-1 proto-oncogene) and several of its known regulated genes were upregulated during injury, possibly reflecting changes in cellular composition. In contrast to widespread changes in gene expression, only a few changes in alternative exon usage could be detected because of injury. There were more than 1,000 changes in retained intron usage, however. Unexpectedly, most of these retained introns have not been described yet but could be validated using direct RNA nanopore sequencing. In addition to changes from injury, our model allowed regional analysis of gene expression. Comparing the segments rostral and caudal to the injury site in naïve animals showed 525 differentially regulated genes and differential regional use of retained introns. We did not detect changes in the serotonin receptor 2C editing that were implicated previously in this spinal cord injury model. Our data suggest that regulation of intron retention of polyadenylated pre-mRNA is an important regulatory mechanism in the spinal cord under both physiological and pathophysiological conditions.

Expanding the Molecular Alphabet of DNA-Based Data Storage Systems with Neural Network Nanopore Readout Processing.

DNA is a promising next-generation data storage medium, but challenges remain with synthesis costs and recording latency. Here, we describe a prototype of a DNA data storage system that uses an extended molecular alphabet combining natural and chemically modified nucleotides. Our results show that MspA nanopores can discriminate different combinations and ordered sequences of natural and chemically modified nucleotides in custom-designed oligomers. We further demonstrate single-molecule sequencing of the extended alphabet using a neural network architecture that classifies raw current signals generated by Oxford Nanopore sequencers with an average accuracy exceeding 60% (39× larger than random guessing). Molecular dynamics simulations show that the majority of modified nucleotides lead to only minor perturbations of the DNA double helix. Overall, the extended molecular alphabet may potentially offer a nearly 2-fold increase in storage density and potentially the same order of reduction in the recording latency, thereby enabling new implementations of molecular recorders.

Genomic regions associated with virulence in Setosphaeria turcica identified by linkage mapping in a biparental population.

Northern corn leaf blight (NCLB) and sorghum leaf blight (SLB) are significant diseases of maize and sorghum, respectively, caused by the filamentous fungus Setosphaeria turcica. Strains of S. turcica are typically host-specific and infect either maize or sorghum. Host specificity in this pathogen is attributed to a single locus for maize and a second distinct locus for sorghum. To identify the genetic basis of host specificity in S. turcica, we generated a biparental population of S. turcica by crossing strains specific to maize and sorghum, phenotyped the population for leaf blight on sorghum and maize, genotyped the population to create a linkage map of S. turcica, and located candidate virulence regions. A total of 190 ascospores from 35 pseudothecia were isolated from the cross of maize and sorghum-specific strains. Greenhouse phenotyping of the biparental population (n = 144) showed independent inheritance of virulence, as indicated by a 1:1:1:1 segregation for virulence to maize, sorghum, both maize and sorghum, and avirulence to both crops. The population and host-specific parent strains were genotyped using genome skim sequencing on an Illumina NovaSeq 6000 platform resulting in over 780 million reads. A total of 32,635 variants including single nucleotide polymorphisms and indels were scored. There was evidence for a large deletion in the sorghum-specific strain of S. turcica. A genetic map consisting of 17 linkage groups spanning 3,069 centimorgans was constructed. Virulence to sorghum and maize mapped on distinct linkage groups with a significant QTL detected for virulence to maize. Furthermore, a single locus each for the in vitro traits hyphal growth rate and conidiation were identified and mapped onto two other linkage groups. In vitro traits did not correlate with in planta virulence complexity, suggesting that virulence on both hosts does not incur a fitness cost. Hyphal growth rate and conidiation were negatively correlated, indicating differences in hyphal growth versus dispersal ability for this pathogen. Identification of genetic regions underlying virulence specificity and saprotrophic growth traits in S. turcica provides a better understanding of the S. turcica- Andropogoneae pathosystem.

Using Genomics to Shape the Definition of the Agglutinin-Like Sequence (ALS) Family in the Saccharomycetales.

The Candida albicans agglutinin-like sequence (ALS) family is studied because of its contribution to cell adhesion, fungal colonization, and polymicrobial biofilm formation. The goal of this work was to derive an accurate census and sequence for ALS genes in pathogenic yeasts and other closely related species, while probing the boundaries of the ALS family within the Order Saccharomycetales. Bioinformatic methods were combined with laboratory experimentation to characterize 47 novel ALS loci from 8 fungal species. AlphaFold predictions suggested the presence of a conserved N-terminal adhesive domain (NT-Als) structure in all Als proteins reported to date, as well as in S. cerevisiae alpha-agglutinin (Sag1). Lodderomyces elongisporus, Meyerozyma guilliermondii, and Scheffersomyces stipitis were notable because each species had genes with C. albicans ALS features, as well as at least one that encoded a Sag1-like protein. Detection of recombination events between the ALS family and gene families encoding other cell-surface proteins such as Iff/Hyr and Flo suggest widespread domain swapping with the potential to create cell-surface diversity among yeast species. Results from the analysis also revealed subtelomeric ALS genes, ALS pseudogenes, and the potential for yeast species to secrete their own soluble adhesion inhibitors. Information presented here supports the inclusion of SAG1 in the ALS family and yields many experimental hypotheses to pursue to further reveal the nature of the ALS family.

Interacting impact of maternal inflammatory response and stress on the amygdala transcriptome of pigs.

Changes at the molecular level capacitate the plasticity displayed by the brain in response to stress stimuli. Weaning stress can trigger molecular changes that influence the physiology of the offspring. Likewise, maternal immune activation (MIA) during gestation has been associated with behavior disorders and molecular changes in the amygdala of the offspring. This study advances the understanding of the effects of pre- and postnatal stressors in amygdala gene networks. The amygdala transcriptome was profiled on female and male pigs that were either exposed to viral-elicited MIA or not and were weaned or nursed. Overall, 111 genes presented interacting or independent effects of weaning, MIA, or sex (FDR-adjusted P-value <0.05). PIGY upstream reading frame and orthodenticle homeobox 2 are genes associated with MIA-related neurological disorders, and presented significant under-expression in weaned relative to nursed pigs exposed to MIA, with a moderate pattern observed in non-MIA pigs. Enriched among the genes presenting highly over- or under-expression profiles were 24 Kyoto Encyclopedia of Genes and Genomes pathways including inflammation, and neurological disorders. Our results indicate that MIA and sex can modulate the effect of weaning stress on the molecular mechanisms in the developing brain. Our findings can help identify molecular targets to ameliorate the effects of pre- and postnatal stressors on behaviors regulated by the amygdala such as aggression and feeding.

Ultra-long DNA molecule isolation from plant nuclei for ultra-long read genome sequencing.

Long and ultra-long read DNA sequencing technologies require high molecular weight DNA with high quality and sufficient quantity, which could be challenging to obtain from recalcitrant plant tissues. We describe a protocol to isolate ultra-long DNA from 12 species for ultra-long read genome sequencing. A suitable nuclei lysis buffer is critical for DNA quality and yield. This protocol will enable individual labs to isolate high molecular weight DNA at a rapid pace with low cost from a variety of plant species. For complete information on the use and execution of this protocol, please refer to: Zhang et al. (2020).

Lasting and Sex-Dependent Impact of Maternal Immune Activation on Molecular Pathways of the Amygdala.

The prolonged and sex-dependent impact of maternal immune activation (MIA) during gestation on the molecular pathways of the amygdala, a brain region that influences social, emotional, and other behaviors, is only partially understood. To address this gap, we investigated the effects of viral-elicited MIA during gestation on the amygdala transcriptome of pigs, a species of high molecular and developmental homology to humans. Gene expression levels were measured using RNA-Seq on the amygdala for 3-week-old female and male offspring from MIA and control groups. Among the 403 genes that exhibited significant MIA effect, a prevalence of differentially expressed genes annotated to the neuroactive ligand-receptor pathway, glutamatergic functions, neuropeptide systems, and cilium morphogenesis were uncovered. Genes in these categories included corticotropin-releasing hormone receptor 2, glutamate metabotropic receptor 4, glycoprotein hormones, alpha polypeptide, parathyroid hormone 1 receptor, vasointestinal peptide receptor 2, neurotensin, proenkephalin, and gastrin-releasing peptide. These categories and genes have been associated with the MIA-related human neurodevelopmental disorders, including schizophrenia and autism spectrum disorders. Gene network reconstruction highlighted differential vulnerability to MIA effects between sexes. Our results advance the understanding necessary for the development of multifactorial therapies targeting immune modulation and neurochemical dysfunction that can ameliorate the effects of MIA on offspring behavior later in life.

Association of Common Genetic Variants in the CPSF7 and SDHAF2 Genes with Canine Idiopathic Pulmonary Fibrosis in the West Highland White Terrier.

Canine idiopathic pulmonary fibrosis (CIPF) is a chronic fibrotic lung disease that is observed at a higher frequency in the West Highland White Terrier dog breed (WHWT) and may have molecular pathological overlap with human lung fibrotic disease. We conducted a genome-wide association study (GWAS) in the WHWT using whole genome sequencing (WGS) to discover genetic variants associated with CIPF. Saliva-derived DNA samples were sequenced using the Riptide DNA library prep kit. After quality controls, 28 affected, 44 unaffected, and 1,843,695 informative single nucleotide polymorphisms (SNPs) were included in the GWAS. Data were analyzed both at the single SNP and gene levels using the GEMMA and GATES methods, respectively. We detected significant signals at the gene level in both the cleavage and polyadenylation specific factor 7 (CPSF7) and succinate dehydrogenase complex assembly factor 2 (SDHAF2) genes (adjusted p = 0.016 and 0.024, respectively), two overlapping genes located on chromosome 18. The top SNP for both genes was rs22669389; however, it did not reach genome-wide significance in the GWAS (adjusted p = 0.078). Our studies provide, for the first time, candidate loci for CIPF in the WHWT. CPSF7 was recently associated with lung adenocarcinoma, further highlighting the potential relevance of our results because IPF and lung cancer share several pathological mechanisms.

The 'in vivo lifestyle' of bile acid 7α-dehydroxylating bacteria: comparative genomics, metatranscriptomic, and bile acid metabolomics analysis of a defined microbial community in gnotobiotic mice.

The formation of secondary bile acids by gut microbes is a current topic of considerable biomedical interest. However, a detailed understanding of the biology of anaerobic bacteria in the genus Clostridium that are capable of generating secondary bile acids is lacking. We therefore sought to determine the transcriptional responses of two prominent secondary bile acid producing bacteria, Clostridium hylemonae and Clostridium hiranonis to bile salts (in vitro) and the cecal environment of gnotobiotic mice. The genomes of C. hylemonae DSM 15053 and C. hiranonis DSM 13275 were closed, and found to encode 3,647 genes (3,584 protein-coding) and 2,363 predicted genes (of which 2,239 are protein-coding), respectively, and 1,035 orthologs were shared between C. hylemonae and C. hiranonis. RNA-Seq analysis was performed in growth medium alone, and in the presence of cholic acid (CA) and deoxycholic acid (DCA). Growth with CA resulted in differential expression (>0.58 log2FC; FDR < 0.05) of 197 genes in C. hiranonis and 118 genes in C. hylemonae. The bile acid-inducible operons (bai) from each organism were highly upregulated in the presence of CA but not DCA. We then colonized germ-free mice with human gut bacterial isolates capable of metabolizing taurine-conjugated bile acids. This consortium included bile salt hydrolase-expressing Bacteroides uniformis ATCC 8492, Bacteroides vulgatus ATCC 8482, Parabacteroides distasonis DSM 20701, as well as taurine-respiring Bilophila wadsworthia DSM 11045, and deoxycholic/lithocholic acid generating Clostridium hylemonae DSM 15053 and Clostridium hiranonis DSM 13275. Butyrate and iso-bile acid-forming Blautia producta ATCC 27340 was also included. The Bacteroidetes made up 84.71% of 16S rDNA cecal reads, B. wadsworthia, constituted 14.7%, and the clostridia made up <.75% of 16S rDNA cecal reads. Bile acid metabolomics of the cecum, serum, and liver indicate that the synthetic community were capable of functional bile salt deconjugation, oxidation/isomerization, and 7α-dehydroxylation of bile acids. Cecal metatranscriptome analysis revealed expression of genes involved in metabolism of taurine-conjugated bile acids. The in vivo transcriptomes of C. hylemonae and C. hiranonis suggest fermentation of simple sugars and utilization of amino acids glycine and proline as electron acceptors. Genes predicted to be involved in trimethylamine (TMA) formation were also expressed.